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Q>^ ■ ABSTRACT 

m • 

fN| Context. One possible way to investigate the nature of the primordial power spectrum fluctuations is by investigating the statistical properties 
of the local maximum in the density fluctuation fields. 

Aims. In this work we present a study of the mean correlation function, and the correlation function for high amplitude fluctuations (peak- 
' peak correlation) in a slighlty non-Gaussian context. 
Methods. From the definition of the correlation excess, we compute the Gaussian two-point correlation function and, using an expansion in 
Generalized Hermite polynomials, we estimate the correlation of high density peaks in a non-Gaussian field with generic distribution and power 
O I spectrum. We also apply the results obtained to a scale-mixed distribution model, which correspond to a nearly Gaussian model. 

Results. The results reveal that, even for a small deviation from Gaussianity, we can expect high density peaks to be much more correlated than 
in a Gaussian field with the same power spectrum. In addition, the calculations reveal how the amplitude of the peaks in the fluctuations field 
^ is related to the existing correlations. 

Conclusions. Our results may be used as an additional tool to investigate the behavior of the N-point correlation function, to understand how 
• • non-Gaussian correlations affect the peak-peak statistics and extract more information about the statistics of the density field. 

• I— I 
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1. Introduction tests of this kind are not a straightforward task, since there is 

no adequate general test for every kind of model. To attack 

Investigation of the statistical properties of cosmological this problem, any effort to better understand how the statistical 

density fluctuations is a very useful tool to understanding properties of the density fluctuation field affect the observed 

the origin of the cosmic structure. Roughly, cosmological Universe is welcome, since it may bring extra pieces of infor- 

models to describe primordial fluctuations can be divided mation to the investigation of cosmic structure, 

in two classes: Gaussian and non-Gaussian. The most ac- Due to the great importance of characterizing 

cepted model for structure formation assumes initial quan- non-Gaussian signatures, many statistical approaches 

tum fluctuations created during inflation and amplified by have been used to study the distribution of fluctua- 

gravitational effects. The standard inflationary models pre- tions in th e cosmic microw av e background radiation 

diets an uncorrected random field, with a scale-invariant (CMBR) (jC hiang et al. 2003'; Bond & Efstathiou 1987; 

power spectrum, which follows a nearly-Gaussian distribution lEriksen et al. 2005.; .Cabella et al. 2004 ; Komatsu et al. 2003) 

(Guth & PI 1982; Gangui et al. 1994). However, non-Gaussian and flie large scale structti re (LSsTtFrith et al. 2005; 2003; 

fluctuations are also allowed in a wide class of alternative mod- IVerde et al. 2 00ll iFrv 1 98.il) . One possible way to investigate 

els, such as: the m ultiple interactive fields (lAUen et al. l"987l the nature of the primordial power spectrum fluctuations is by 

ISalopek et al. 19891) . the c osmic defects models ( Kibble 1976) investigating the statistical properties of the local maximum 

and the hybrid m odels ( Magueiio & Brandenbergher 2000i: in the density fluctuation fields. Since some of the peak 

iBattve & Weller 2000) . By discriminating between different properties, such as number, frequency, coiTelation, height 

classes of models, the statistical properties of the fluctuations and extrema, are highly dependent upon the statistics of 

field can be used to investigate the nature of cosmic structure, the fluctuation field, we can gather information about the 

However, non-Gaussian models comprehend an infinite range statistical distribution function by studyi ng the morphological 

of possible statistics. As a consequence, performing statistical properties of the fluctuations fields (iBardeen et al. 198^ 



2 



A.P.A. Andrade et al.: High Order Correction Terms for The Peak-Peak Correlation Function in Nearly-Gaussian Models 



iNovikov et al. 1999t iLarson & Wandelt 20051) . Other use- 
ful statistical estimator applied to investigate density 
fluctuations are the wavelets tools ("McE wen et al. 2004: 
^ukherjee & Wang 2004; Vielva et al. 200 4J, the phase 
correlations dColes et al. 20041 iNaselskv et al. 20051) and 
the most widely used estimator: the N-point correla- 
tions in phase (L ewin et al. 1999; Verde & Heavens 2001; 



Kgmgtsu^^^lj^^QS^) or density spaces ( Takad^j^J^iii2223 

Coles & Jones 199l1 iHeavens & Sketh 19991 

Heavens & Gupta 200ll lPeebles 200lH Eriksen et al. 200^! 

The extensive use of the two-point correlation function to 
characterize the statistical properties of the fluctuations field 
is justified by the mathematical simplicity in the Gaussian 
condition, since the mean correlation function can be ob- 
tained analytically and it completely specifies a Gaussian 
distribution (as weU as its power spectrum). However, this 
assumption is not true for a non-Gaussian case, where higher 
order correlations may give a significant contribution, de- 
spite the great effort demanded to compute a wide range 
of correlations. However, it is possible to detect primordial 
non-Gaussianity with a non-zero measure of the N-point 
correlation (n > 3). In order to achieve a good descrip- 
tion of non-Gaussian signatures in cosmic structure, many 
works have been done to estimate the three-point correla- 
tion function and the related bi or tri-spectrum for a few 



Komatsu et al. 2003; 
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Eriksen et al. 2005; Komatsu et al. 2005 


Gaztanaga & Wagg 20Q4|; 


(Hansen et al. 2004 



that, with the advent of the CMBR experiments and the high 
quality surveys in cosmology, the N-point correlation function 
will be the main statistical descriptor for the cosmic structure. 

In this work we present a technique to extract non-Gaussian 
components from a two-point correlation function. In presence 
of a non-Gaussian component, we point out that even between 
two points it is possible to estimate the influence of higher 
order correlations for models with different statistical descrip- 
tion. By performing the calculation of the high order correction 
terms for the peak-peak correlation function, we show how the 
amplitude of peaks is dependent of the correlations involved. 
We point out that for a non-Gaussian statistics the higher order 
correlations impose some restrictions to the amplitude of high 
density peaks. 

This paper is divided as follows: in Section 2 we give a 
general description of a random variable field and estimate the 
peak-peak correlation function for a Gaussian field. In Section 
3 we describe the general treatment to obtain the peak-peak 
non-Gaussian correlation function. In Section 4 we apply the 
calculations for a slightly non-Gaussian model in two steps: 
first we consider an approximated solution for a non-Gaussian 
model with null high order correlations and finally we estimate 
the complete solution for a slightly non-Gaussian field. In sec- 
tion 5, we summarise our results and discuss the possibilities to 
use the peak-peak correlation function as a statistical descriptor 
of the density fluctuation field. 



2. Random Variable Fields 

Most of the models for the early universe (i.e. inflation) actually 
predicts the fluctuation field to be random. This requires that 
6(x) can be treated as a random variable in the 3D-space and 
the assumption that the universe is a random realization from a 
statistical ensemble of possible universes. 

We define a random variable, 6, using the fact that, instead 
of knowing its exact value, we only know how to measure var- 
ious values of 6i, 62, ... 6„, which define a random variable 
field, under certain experimental conditions. Therefore, a ran- 
dom variable can only be characterized by a certain statisti- 
cal ensemble of realizations. When we say that a random vari- 
able is known, it means that we only know the statistical sam- 
ple which characterizes it. To completely describe the statis- 
tical properties of a random variable, S, we define the proba- 
bility density function, P[6], which can be obtained from the 
Fouri er transform of the characte ristic distribution function, 0^ 
dGnedenko & Kolmogorov I968I) : 



1. r 



(1) 



The characteristic function, ©^(m), can be obtained by the 
McLaurin series for the moments, m: 



v-i (iu)" 

&s(u) = 1 H- > — — OT„, for : nin 



(2) 



Another possibility to obtain the characteristic function of a 
random variable is to use the distribution in series of cumulants: 



0^(m) = exp 



n=\ 



for : kn = {{6- < 6 >)"} 



(3) 



Since the physical importance of the cumulants k„ in Eq.|3lde- 
creases as n increases, it is usual to confine the statistical cal- 
culations of random variables to the first few terms of the cu- 
mulant distribution series and, for convenience, set the higher 
order terms to zero. However, the calculations presented in the 
next sections show that, even if the cumulants terms are very 
small, they can significantly contribute to the statistical descrip- 
tion of non-Gaussian fields. 

2.1. Correlations in a Random Field 

The main numerical indication of the correlation degree be- 
tween random variables are the N-order correlation functions. 
The autocorrelation (or double correlation) for a random vari- 
able 6 is defined by: 



(4) 



The triple correlation is similarly defined in terms of all possi- 
ble combinations between the three variables, being: 



K3[6id2S3] = (61626^) - {Si}K[6263] - {62} K[6i6i] 

-{6i}K[6iS2]-{6l}{62){63} (5) 
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Note that, in the case where we have the same three vari- 
ables (^1 -62 - Si), the three-point correlation function is sim- 
ilar to the third cumulant of the distribution. Higher order cor- 
relations between several variables can be defined, in a simi- 
lar way, by the difference between all possible correlations in- 
volved. 

For a statistical process where the correlation functions of 
order greater than one are null, we set the variable described by 
this correlation function as not random, or deterministic. In the 
case where the correlation functions of order greater than two 
are null, we have a Gaussian variable. For the case of corre- 
lation functions of order greater than two not completely null, 
the variable is considered to be non-Gaussian. In this sense, we 
can say that a Gaussian random field is a simplified version of 
a general random field. 

Usually, the cosmological density fluctuation field is sta- 
tistically described by the mean correlation function, ^(r), ap- 
plied to a galaxy or a cluster distribution, with two-point mean 
separation defined by r (Peebles 1980). For an isotropic and 
homogenous field, the correlation function is defined as the ex- 
cess of probability for a density field described by a Poisson 
distribution. Therefore, the probability to find two points, in a 
volume dV]dV2, separated by a distance ri2 is given by: 



dP = n^dVxdVi [1 +^(^12)] 



(6) 



Describing the fluctuations field in Fourier modes, we have: 



which is equivalent, in a continuous space, to: 

_V 

(2^ 



(7) 



(8) 



2.2. High Density Peaks in a Gaussian Random Field 

For a Gaussian random field, the n-dimensional probability 
density function can be estimated from the Fourier transform 
of the characteristic function (Eq.^andlJJ defined for moments 
distribution of s <2. 



P'^[6u62-6„] ^ dMiJo dM2-Jo du„ 



exp 



(Ittt Jo 

2 n 
^ — ^ KASa, 5l])UaUp 



s=\ ' afi=\ 



exp 



(9) 



The expression above can be reduced to: 



For a bidimensional Gaussian fluctuation field with zero 
mean, the correlation matrix can be obtained and inverted, re- 
sulting in: 



1 I 1 V- 



e 



(11) 



where A = , for: 

& = {Si52), r^\x2- xi\, o-^ = (52) = (52), and: 



(12) 



where is equal to ^^2(^1 , ^2), the mean correlation function 
for a Gaussian field. 

To find the correlation function between high density peaks, 
we calculate the probability of i5i and 62 to exhibit density 
values which are larger than the variance field cr, by a factor 
rjirj > 0) . For a Gaussian fiel d, this probability is given by the 
integral jPadmanabhan 199^ : 

P^[6i>j]o-,62>w]^ I ddi I d62P^[6i,S2]- (13) 

The integral above can be obtained by substituting Eq. |5] 
into Eq.^l For a weakly correlated field, where A - ^ « 1, 
and high density peaks, 77 > 1, the integral above will be: 



P^[6i > Tjo-,62 > w] 



f d6i r d62 

Jno- Jno- 



p=0 



(2A)" 



where F is the Err function (erf): 
27r 

The final result is 



APer'^''''\dA. 



P [61 > ria-,62 > w] - 



1 1 



(1 +2Ar]^). 



(14) 



(15) 



(16) 



Redefining expression the probability to find peaks 61 
and 62 with density rj times the variance, cr, is: 



(17) 



P''[6u62...6n] = -\^\det\\K2(6a,6My^- 

[In) 2 L J 

n 

^ [6a - Kl{6a)] [dfS - Kiidp)] 



exp 



(10) 



where det\\K2{5a,Sp)\\ and Wa^pW are, respectively, the determi- 
nant and the inverse of the correlation matrix. 



where [P [S > rjcr)]] is the mean probability to find high den- 
sity peaks in a bi-dimensional random field, and [1 + ^{r)} is 
the probability excess expressed in terms of the mean correla- 
tion function, Then, for weakly correlated fields: 



e^{r) ^ lAif ^ e^^'f -l^exp 2 (^)' ^, 



- 1. 



(18) 
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3. Correlations in a Generic Non-Gaussian Field 

The coiTelation function of high density peaks in a non- 
Gaussian field can also be computed using Eq. [O] except that 
we have to consider the non-Gaussian probability of finding 
both peaks di and 62, the P'^'-' [6 1,62]- This probability can be 
estimated using Eq. |9l performing the sum over terms of or- 
der higher than 2. Since the summation is also computed for 
additional terms, the result can be synthesized by; 



which results in: 



(19) 



[Si > rio-,62 > rjcr] 
1 1 



(20) 



where the factor C^*^ carries the corrections terms of higher 
order (s > 3). 

One intriguing question we could ask is how important the 
high order residual terms for a slightly non-Gaussian statistics 
are . First, we could consider the approximated case in which 
the extra calculation in Eq.|5]were avoided, so correlation terms 
of order (s > 2) could be neglected {C^'^ ~ 0) and the non- 
Gaussian peak-peak correlation function would be reduced to: 



(r) ^ exp 



1, 



(21) 



This solution is similar to the calculation for a Gaussian 
random field, except for the fact that we were considering a 
non-Gaussian mean correlation function. 

At this point we want to assess the robustness of the as- 
sumption C^'^ ~ 0, stated in the previous paragraph. For this 
purpose, we compute the higher order correction terms consid- 
ering a slightly non-Gaussian component, this means a very 
small contribution to correlations of order greater than two. 
However, the non-Gaussian probability described in Eq. |9]can 
not be calculated using the Fourier transform of the character- 
istic function, since there is no analytical solution for that ex- 
pression. One possible way to obtain P'^'^[Si,62], as suggested 
by (Gnedenko & Kolmo2orov 1968), is to expand it in a series 
of Generalized Hermite polynomials, H: 



P'''^[6u62...6„]^P''[6u62...6„] 

1 + ^ — ^ bs(6a, 6p...6JHa^p,„^ [6 - Ki(6)] 



5' 

-3 ' a£...(D=l 



(22) 



where bs is the quasi-moment function and ki is the first cumu- 
lant of the distribution. The definitions of bs and H are given in 
Appendix A. 

By Eq. |52 the non-Gaussian probability is expressed in 
terms of a Gaussian probability added to higher order {s > 



3) correction terms, which are related to the deviation from 
Gaussianity. Combining equations ^] and |22lwe can es- 
timate the higher order corrections term for a bi-dimensional 
using the following: 



CO ^ 2 



s=3 a^...oj=l 



(23) 



Note that the calculation of C^"^ starts at s > 3, ensuring 
that the terms used in the expansion are related to the non- 
Gaussian contributions. Details of the calculation involved in 
this higher order terms are also presented in Appendix A. 

4. Correlations in a Nearly-Gaussian Field 

4.1. The Mixture Model 

In order to estimate how the high order terms affect the corre- 
lation between high density peaks, we estimated the for a 
Gaussian and a slightly non-Gaussian field, computing the ap- 
proximated solutions (Eq. I21li and the full calculation of the 
expansion (Eq. I22> until s < 6 order. For this comparison, we 
have considered a mixed probability distribution, as proposed 
by Ribeiro, Wuensche & Leteher (2001), hereafter (RWL). 

The general procedure to create a wide class of non- 
Gaussian models is to admit the existence of an operator which 
transforms Gaussianity into non-Gaussianity according to a 
specific rule. An alternative approach to study non-Gaussian 
fields was proposed by RWL, in which the PDF is treated as 
a mixture: P{iJ/) - (1 - a)fi((p) + af2{(p), where f]{<p) is a 
(dominant) Gaussian PDF and /2(0) is a second distribution, 
with < a < 1 . The a parameter gives the absolute level of 
Gaussian deviation, while /2(0) modulates the shape of the re- 
sultant non-Gaussian distribution. RWL used this mixed sce- 
nario to probe the evolution of galaxy cluster abundance in the 
universe and found that even at a small level of non-Gaussianity 
(a ~ lO"'* - 10"^) the mixed field can introduce significant 
changes in the cluster abundance rate. 

The effects of such mixed models in the CMBR power 
spectrum, combining a Gaussian adiabatic field with a sec- 
ond, non-Gaussian isocurvature fluctuation field, to produce 
a positive skewness density field was discussed by Andrade, 
Wuensche & Ribeiro ( |2004 (hereafter AWR04). In this ap- 
proach, they adopted a scale dependent mixture parameter and 
a power-law initial spectrum to simulate the CMBR tempera- 
ture and polarization power spectra for a flat A-CDM model, 
generating a large grid of cosmological parameters combina- 
tion. The choice of a scale-dependent mixture is not unjusti- 
fied, since it could fit both CMB and high-z galaxy clustering 
in the Universe (e.g. AWR04; Avehno & Liddle 2004; Mathis, 
Diego & Silk 2004). At the same time, in mixed scenario, the 
scale-dependence acts in order to keep a continuous mixed 
field inside the Hubble horizon. Simulation results show how 
the shape and amplitude of the fluctuations in CMBR depend 
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upon such mixed fields and how a standard adiabatic Gaussian 
field can be distinguished from a mixed non-Gaussian one. 
They also allow one to quantify the contribution of the sec- 
ond component. By applying a test on recent CMBR ob- 
servations, the contribution of the isocurvature field was es- 
timated by Andrade, Wuensche & Ribeiro (2005) (hereafter 
AWR05) as a^^ 0.0004 ± 0.00030 with 68% confidence Hmit. 
In the present work, we also investigate the predictions of scale- 
dependent mixed non-Gaussian cosmological density fields for 
the peak-peak correlation function. 

4.2. The Two-Point Correlation Function 

To obtain the mean correlation function, ^r, we have computed 
the Fourier transform of the power spectrum related to a pure 
Gaussian field and a mixed non-Gaussian PDF. In this sense, 
we rewrite Eq.|S] which is equivalent to: 



V 



(27r)3 



P(k)k^dk 



J, "^Jv 



-ikr 



(24) 



10 
10' 



10 

10'= 







- Gauss 

- Log-Normal; oq^IxIO""' 
Mix(G+LN); ao=lxlO"^ 
Mix(G+LN); ao=3xlO"3 


'■ \ 
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Fig. 1. Tlie mean two-point correlation function estimated for non- 
correlated or weakly correlated fields: pure Gaussian, pure log-normal 
and Mixed field with ao ~ 10"^. 



being u - cosO. For an isotropic field, we have: 



" (2^ Jo 



kr 



k^dk 



(25) 



For the mixture correlation function, we consider a mixed- 
scale power spectrum described as: 



P'^'-'(k) = A„P(k)[l + M(aQ)k] = /3k"[l + M(aQ)k] 



(26) 



This power spectrum was estimated by the correlated 
mixed-model which considers a possible mixture, inside the 
horizon, between fluctuations of an adiabatic Gaussian field 
and an isocurvature non-Gaussian one (AWR04). In this model, 
ao is the mixture parameter that modulates the contribution of 
the isocurvature field (ao - (0.00042 ± 0.00030)) as estimated 
from recent CMB observations. For a null ao, the field is purely 
Gaussian with a simple power law spectrum. A„ is the normal- 
ization constant of the primordial power spectrum, estimated as 
~ 1.3 X 10"" for recent CMBR observations (AWR05). M(ao) 
is the mixed term, which accounts for the statistical effects of 
the second component in the power spectrum, expressed as a 
functional of both distributions, fj and /2: 



M"'"(ao 



,k) = j[(l- 



aok)fl{v) + aokf2(v)]v dv 



(27) 



Evaluating the integral in Eq.|25lfor a mixed power spec- 
trum, we find an expression for the mean correlation function: 



^^"(/?)=A„ 



1.76 X 10-2 ^ 5.76 X lO-'^M(ao) 



3- (I) 



(28) 



where Rq is about 25 Mpc, the mean correlation width for 
galaxy clusters. 

In Figure 1, we show the mean correlation function esti- 
mated for a pure Gaussian field, M(ao) - 0, for a mixed PDF, 
M(ao) = "^0 (y - l)j and also the individual contribution of 



the non-Gaussian field for ao ~ 10^^. In this plot, it is possi- 
ble to observe the importance of even a small contribution of 
a the second component to the mean correlation function. For 
I- < 1 X 10"^ the non-Gaussian component dominates the cor- 
relation function. This behaviour illustrates the excess of power 
in small scales, as observed in the CMB angular power spec- 
trum in the mixed context (AWR04). 

Inserting the mean correlation function described in Eg. 1281 
into the approximated expression of the correlation function 
for high density peaks (Eq. 12 U for a few classes of PDFs, we 
estimate the functions plotted in Figure 2. In this plot, we 
show that the effect of the second component is still observed, 
and that the correlation function for high density peaks is also 
sensible for different non-Gaussian distributions. Comparing 
Figures 1 and 2, we see how the high density peaks can be 
much more correlated than the mean field for a non-Gaussian 
case, especially in small scales, ^ ~ 10"^ where is nearly 
two orders of magnitude greater than the mean correlation. 

With the help of a program that performs algebraic and nu- 
merical calculations, we actually computed C^'^, as indicated 
in Eq.|53 for correlations up to the 6-th order (.? < 6). This limit 
was set in order to keep a meaningful non-Gaussian distortion, 
avoiding more time consuming calculations. We do not present 
in this section the full computed expression since it contains 
hundreds of non-linear terms in the mean correlation function 
i^r) and 77, where the coefficients are the high order correlations 
(or cumulants). In fact, the C^'^ is more accurately described 
as C^'^i^r)- In Appendix A, we show the steps to calculate 
the quasi-moment function, bi,„, and the Hermite Polynomials, 

Him- 

In general, correlations of y erv high order tend to zero 
dCnedenko & Kolmo^orov 1968h . the most extreme case being 
the normal distribution, where all cumulants of > 3 are null. 
Deviations from Gaussianity are set by the increment of non- 
vanishing cumulants in the expansion of C^'~' ■ A possible ques- 
tion that may be raised is related to the convergence and nor- 
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Fig. 2. The two-point correlation function computed for 2cr density 
peaks in both a pure Gaussian and mixed context (q-q ~ 10"^). For 
this estimation, we used: M'^""'(ao) = 0; M^^'^Co-o) = Oq; M'^'"'(ao) = 
ffo(3Vf - i);M^''(q'o) = 2ao andM^'^(Q'o) = ao{^- l). 



malization of the expansion in Eq.|23lin which some terms are 
set as non-vanishing. However, investigation of general nearly 
gaussian deviations have already been performed by the use of 
the Edgeworth expansion in one-dimension by setting the func- 
tion to z ero in higher orders terms and normalizing it appro- 
priately ( Martinez-Go nzalez et al. 20021) . Following these au- 
thors, we have considered, as a work hypothesis, the following 
behavior for such coefficients: we have set the terms .? > 7 to 
zero and the cumulants involved in the quasi moment function 
to 10"^ for 3 < i < 6. 

A very interesting result obtained is summarized in Figure 
3. For this plot, in order to follow the absolute behavior of the 
corrections terms in the nearly Gaussian field, we have set the 
normalization value of the mean correlation function to one, 
(^r ~ 1) and explore the rj dependence of C^'^i^r) with the 
purpose to understand how the amplitude of the fluctuations is 
related with the high order correlations. While computing only 
third-order terms, the describes a gradual enhancement in 
the correlation function. For fourth and sixth-order terms, the 

and C^, in the case of small Gaussian deviation, what is 
observed is that the CjJ and just overlap each other, mean- 
ing that the fifth order correlations do not contribute signifi- 
cantly. However, the fourth order correlation describes peaks 
of maximum density allowed by a weakly correlated field at 
about J] ~ 1.8 cr. For higher order correlations, the Cf^ shows 
a very large increment in the two point correlation for densi- 
ties about 2.5 cr, and a nearly null contribution to peaks with 
amplitudes {rj > 3cr). 

In this sense, we conclude that one should not expect very 
high density peaks for the specified, slightly non-Gaussian, cor- 
related field. However, factors or order 3 < .? < 6, are quite 
significant for peaks with amplitude up to 2.5 cr and are too 
far from a null correction, even if we consider weakly corre- 
lated fields. Then, we can not consider the simplified solution 
in Eq. |^ as a good description of the two-point correlation 



function for high density peaks in despite of our choice of con- 
sidering just a small deviation from Gaussianity. It is important 
to note that this result is independent of the power spectrum or 
the mean correlation function. It only shows the influence of 
higher order correlations in the amplitude of the field of fluctu- 
ations. 

Analyzing the relation between C,j and rj we gain some in- 
sight about the amplitude of such high order terms, since C,, 
controls the amplitude of the permitted peaks in the fluctuations 
field. Increasing values of congelations with s > 6 imply a high 
probability of very high amplitude (very rare) peaks, which 
contradict the observations of large scale structures. However, 
when we impose correlation levels of the order 10"^' up to sixth 
order, we favour the existence of peaks up to 3cr, what is very 
reasonable for a nearly Gaussian field. 

In Figure 4 we show the behaviour of the two-point correla- 
tion function estimated for a Gaussian, for a mixed ap- 
proximated solution, and for the non-Gaussian complete 
solution, ^^'-^ + C^'^i^r) ■ In this plot, we have set the amplitude 
threshold of 77 ~ 2 and test the dependence on ^r- The observed 
eff'ect of C^'^i^r) is to amplify the correlations between high 
density peaks. This result is valid for the case of an increasing 
mean correlation function, and has no dependence of the mixed 
model. While non-Gaussian deviation tends to add non-vanish 
high order correlations, we conclude that we can expect high 
density peaks to be much more correlated even in a slightly 
non-Gaussian model. 

5. Discussion 

In this work, we have estimated the two-point mean coiTelation 
function and the peak-peak correlation function for the den- 
sity field. In the Gaussian case, the calculations are simplified, 
since the Fourier transform of the power spectrum completely 
describes the random variable. However, for a non-Gaussian 
field, the calculations are much more complicated, since high 
order correlations between these two-points may not vanish and 
strongly contribute to the final function, even in a small devia- 
tion from the Gaussian case. 

In this work, we showed that, when considering a mixed 
model, both the mean coiTelation and the peak-peak correla- 
tion functions are much more intense in small scales than in 
the Gaussian case. This result can be particularly relevant since 
it is generally accepted that galaxies form in high density re- 
gions. In addition, we conclude that the peak-peak coiTelation 
function is quite sensitive to the PDF of the fluctuations field, 
especially for a mixed model. This result suggests that it is pos- 
sible to use the peak-peak correlation function as a test for the 
nature of cosmic structures. Nevertheless, we have to be care- 
ful when approximating terms for high order correlations, since 
the peak-peak correlation function is very sensitive to the cor- 
rection terms. We also point out that correlations of order > 2 
can be a very important tool to characterize non-Gaussian fields 
and definitely deserve a deeper investigation. Estimating high 
order correlations allows us to investigate the behaviour of the 
N-point correlation function, as well as gather more informa- 
tion about the amplitude of the expected high density fluctu- 
ation field. As seen in Figure 3, correlations may restrict the 
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Fig. 3. The behaviour of the higher order factor, C'^'^, for a fixed mean 
correlation (^,- ~ 1). The curves in A) show the factor estimated for 
third order correlation, C^, fourth, CjJ, and fifth, C,^. In B) the sixth 
order factor, C*, is also plotted. 



Fig. 4. The behaviour of the two-point correlation function is shown 
for three estimated cases: a pure Gaussian, f^"'"' (lower curve, or- 
ange), a mixed approximated solution, with ao ~ 3 x 10"^ 
(mid curve, black); and for the complete non-Gaussian solution, 
+ C^°(^r), estimated for a PDF of the type: (Gauss + Exp) with 
tto ~ 3 X 10^' 



understanding, the investigation of the peaks statistics in CMB 
datasets should be realized over several datasets and different 
frequencies. Indeed, one possible manner to minimize the fore- 
grounds effects is to analyse the most sensitive and cleaned 
map, such the WMAP three-year di -biased internal linear com- 
bination, WMAP-DILC map, or the WMAP coadded map, 
with the combination of Q-i-V-i-W frequencies. This analysis 
will be the next test applied to the non-Gaussian mixed model. 
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amplitude of the density peaks. Furthermore, from the ampli- 
tude of the peaks found in the density field (in CMB or LSS 
fluctuations) we are able to extract more information about the 
statistics of the density field. It is also good to remember that 
the presence of correlations of order s > 3 lead to formation of 
structures in earlier times than would be expected for a model 
with the same power spectrum but with weaker spatial correla- 
tions. The resuls presented in this paper may be used to set new 
constraints in structure formation models. 

One possible application of this method on investiga- 
tions of primordial non-Gaussianity could be implemented 
in the search for maximum amplitude fluctuations of the 
full sky CMB temperature maps, such those derived from 
WMAP (^Spergel et al. 1989) and the future Planck mission 
(I Wright E.L. 200Q) . However, variations in the number of den- 
sity peaks and their correlation may as well be related to non- 
Gaussian Galactic foregrounds or other contaminants. Once 
detected such non-Gaussian trace, we have to be very care- 
ful before assigning it to a primordial origin. To avoid mis- 
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sible way to obtain the quasi-moment functions, bs, is relating 
to the correlation function, since: 



exp 



CO n 



2 K,{Xa,Xi3,...,X^)Z 

.9=3 a,p,...,cij=\ 



(A.l) 



.v=3 (2 ,/?,..., tL»=l 



The Generalized Hermite polynomials can be obtained by the 
definition: 



HaB...<Ax\ = e' 



,0W 



dxajx dXf 



dXa 



(A.2) 



where: (f>[x] - \YI'ap=\ dapXaXp, being \\aap\\ the elements of 
the inverse correlation matrix ||A;2(fa, f/3)||. The unidimensional 
case of Eq. |23is known as Edgeworth series, while the bidi- 
mensional case corresponds to: 



(A.3) 



1 + ^ ^ ^ bimHi„, [6 - Ki (6)] 

s=3 l+m=s 

where: bi,„ = bi+„,[Si-^, 62-^62], 
I times m times 
1;2_— 2- 
/ times „, times 



Performing the calculations for the quasi-moment function 
we have: 

bi^K[di,62,63]; 
b4^K[6u 62,63,64]; 
bs^K[6i,d2,6i,64,65]; 

be = K[6i,62, 6i, 64, 65, de] + 10 [K3[6i,62, d^Wdd^, S5, Se]] . 

Note that the first five terms of bs are equivalent to the 
correlation function, only the sixth term have additional terms. 
As an example, for 5 = 3, we obtain four different forms of by. 



b3Q^K3[6i,6uSi,]^{6])- 
bo3=K3[62,62,62,] = {sl); 
b2i =K3[6u6uS2,] = {siS2); 
bn = K3[6u62,62,] = {Si6l). 



Appendix A: Calculation of and H 

To obtain the non-Gaussian probability for a multi-dimensional 
case we have to perform the calculation in Ea.l22lexpanding the 
Hermite Polynomials and the quasi-moment function. One pos- 



Proceeding in a similar manner, we find five terms for Z74, 
six for bs and seven terms for bd- 

The best way to find the Hermite Poly nomials is using ex- 
pression dGnedenko & Kolmogorov 1968 



ya = ya[x] 



d(p[x] 

dXa 



lapxp, being : 



(A.4) 
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Ha = Ja, 

Hap — Jo-J/S ~ l^typj 

HaPy - >'ff>'/S>V ~ <^apyy " l^ayj^ — (lypja 

= yaypyy - [aapyy]^ , 

HaPyS = yaypyyyd - \aapyyys\^ - \aapays\^ , 

HaPySto = yaypyyysy<j - {aapyyysycJf^Q - {aapaysyoj}^^, 

Hapy6cM = yaypyyysy,0e - \aapyyysyuiye\^^ - 

\ciap(ly5y(i)ye\^^ - \aap(lyd(l^^^ . 

where we have used the notation {j'jj to indicate a simetrization 
set. 

To perform the calculation above for bim and Him up to (/ + 
m = 6) and perform the integration in d6\ and d62, we used a 
software for algebraic and numerical calculations. 



